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Elastic scattering cross sections of the 89 Y(a,a) 89 Y reaction have been measured at energies E c . m . 
= 15.51 and 18.63 MeV. The high precision data for the semi-magic N — 50 nucleus 89 Y are used to 
derive a local potential and to evaluate the predictions of global and regional a-nucleus potentials. 
The variation of the elastic alpha scattering cross sections along the N = 50 isotonic chain is 
investigated by a study of the ratios of angular distributions for 89 Y(a,a) 89 Y and 92 Mo(a,a) 92 Mo 
at E c . m . ~ 15.51 and 18.63 MeV. This ratio is a very sensitive probe at energies close to the Coulomb 
barrier, where scattering data alone is usually not enough to characterize the different potentials. 
Furthermore, a-cluster states in 93 Nb = 89 Y ® a are investigated. 
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I. INTRODUCTION 

Alpha-nucleus potentials are basic ingredients for the 
calculation of reaction cross sections with a particles in 
the entrance or exit channel. These reaction cross sec- 
tions are included in the calculation of stellar reaction 
rates in nuclear astrophysics which have to be determined 
at stellar temperatures corresponding typically to sub- 
Coulomb energies for reactions involving a particles. 

In several astrophysical applications — such as mod- 
eling the nucleosynthesis in explosive scenarios like p 
process — the reaction rates are taken from statistical 
model calculations [3, HI- These calculations utilize global 
alpha-nucleus optical potential parameter sets. Consider- 
able efforts have been devoted in recent years to improve 
the knowledge of the alpha-nucleus optical potential 0- 
6] . The extensive use of the statistical model calculations 
requires further experimental tests for the global param- 
eterizations. 

The optical potential combines a Coulomb term with 
a complex nuclear potential, composed of real and imag- 
inary parts. The variation of the potential parameters of 
the real part as a function of mass and energy is smooth 
and relatively well understood 0]. On the contrary, the 
imaginary part of the optical potential is strongly energy- 
dependent especially at energies around the Coulomb 
barrier. Therefore, tests of global a-nucleus potentials 
need to focus on experimental information at energies as 
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close as possible to the astrophysically relevant energy re- 
gion. One possible way of testing the different potential 
parameterizations is to carry out alpha elastic scattering 
experiments and compare the measured angular distri- 
butions to the corresponding predictions from the global 
potential parameterizations. However, at astrophysical 
energies the alpha-nucleus elastic scattering cross section 
is non-diffractive and dominated by the Rutherford com- 
ponent. Therefore, the experiments have to be carried 
out at slightly higher energies with high precision. From 
the analysis of the measured angular distributions the 
parameters of the potential can be derived and have to 
be extrapolated down to the astrophysically relevant en- 
ergy region where the relevant alpha-particle induced- 
reactions are taking place. 

The present study focuses on the determination of op- 
tical potentials from elastic scattering cross sections. It 
is an interesting and still open question whether reaction 
cross sections can be predicted from such a potential that 
has only been adjusted to scattering data. But it remains 
as a final goal to find a perfect potential that is able to 
reproduce elastic scattering and reaction cross sections 
and furthermore properties of cluster states like excita- 
tion energies and decay widths. 

Several alpha elastic scattering experiments on even- 
even nuclei 92 Mo, 106 Cd, 112 > 124 Sn, and 144 Sm have been 
performed at ATOMKI in recent years ^4il|. in order 
to investigate the behavior of alpha-nucleus optical po- 
tentials. This work presents the elastic scattering exper- 
iment performed on the 89 Y nucleus to study further the 
systematic behavior of the optical potentials at low ener- 
gies. In all of these cases complete angular distributions 
have been measured at energies close to the Coulomb 
barrier. The chosen energies were low enough to be close 
to the region of astrophysical interest and high enough 
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that the scattering cross section differs sufficiently from 
the Rutherford one. 

Previous studies have focused on semi-magic even-even 
nuclei with N = 82 ( 144 Sm), N = 50 ( 92 Mo), and Z = 50 
( 112,124 Sn). As a natural extension of previous work, 
we are presently working on nuclei that arc cither non- 
magic or even-odd nuclei. Thus, a main motivation for 
our latest experiment on the non-magic nucleus 106 Cd 
was to analyze the influence of shell closures on the a- 
nucleus potential . As a continuation of the systematic 
study of the behavior of the optical potentials, this work 
presents the elastic scattering experiment performed on 
the neutron-magic (N = 50), proton-odd (Z = 39) nu- 
cleus 89 Y. Angular distributions have been measured at 
E c .m.= 15.51 and 18.63 MeV, just above the Coulomb 
barrier (the height of the Coulomb barrier for the 89 Y 
® a system is about 15 MeV). At these energies a reli- 
able test for the global parametrization is possible using 
the new high precision data. Furthermore, we can use 
data from literature [IHIIH to investigate the variation of 
the imaginary part of the optical potential between E c m . 
ps 15.5 and 23.9 MeV in small steps of approximately 2 
MeV. In addition, experimental ang ular distributions are 
available at higher energies [144l7j . A local optical po- 
tential for 89 Y can be established from these data over a 
wide range of energies. This was a further motivation to 
study a scattering on 89 Y. 

A global alpha-nucleus optical potential must be able 
to provide a correct prediction for the elastic scatter- 
ing cross section, and to describe the variation of the 
angular distributions along isotopic and isotonic chains. 
This is especially important for the extrapolation to 
unstable nuclei where the potential cannot be derived 
from experimental scattering data. Galaviz et al. [l(| 
measured the elastic scattering cross sections of the 
112 - 124 Sn(a,a) 112 ' 124 Sn reactions at E Q = 19.5 MeV over 
a broad angular range with small uncertainties of about 
ps 3-4%. The study of both the proton- and neutron- 
rich stable tin isotopes provided important information 
about the changes of the potential parameters with the 
neutron number along the Z = 50 isotopic chain. The 
ratio of the measured cross sections showed an oscillation 
pattern which was very pronounced in the backward an- 
gle region. The analysis of this ratio provides a further 
constraint for global a-nucleus potentials, in particular 
for the extrapolation to unstable nuclei. It needs to be 
highlighted that all global a-nucleus optical potentials 
failed to reproduce either the amplitude or the phase of 
this oscillation pattern in the ratio of the measured angu- 
lar distributions of 112 Sn and 124 Sn at backward angles 
[l(|. This oscillation feature of the experimental data 
was first observed using our precise data on 112 Sn and 
124 Sn, but it was hidden previously because the typical 
uncertainties of the elastic scattering cross sections were 
of the order of 10-15% 0. 

As an extension of our previous work on the isotopic 
chain at Z = 50, now we investigate the behavior of 
the optical potential parameters along the TV = 50 iso- 



tonic chain. We have measured a scattering data of 89 Y 
with small uncertainties over the whole angular range. 
This is a prerequisite to study the ratio of the Ruther- 
ford normalized cross sections of 89 Y(a,a) 89 Y from this 
work and 92 Mo(a,a) 92 Mo previously studied by Fiilop et 
al. H| at Ec. m . = 13.20, 15.69 and 18.62 MeV. Further 
studies could be done on the N = 50 nuclei 86 Kr, 87 Rb, 
and 88 Sr. Low-energy data for the N = 50 nucleus 90 Zr 
at 15 MeV [l9| were already analyzed in 0|. 

This paper is organized as follows. In Sec. HI] we de- 
scribe our experimental procedure. In a first analysis a 
local optical potential is derived from the new experimen- 
tal data and the available scattering data in literature, 
and furthermore bound state properties of cluster states 
in the nucleus 93 Nb = 89 Y <E> a are studied (Sec.|niJ). The 
measured angular distributions as well as the ratio of the 
elastic scattering cross sections of the two N = 50 nuclei 
89 Y and 92 Mo are compared to predictions using sev- 
eral global optical potential parameterizations in Sec. lIVI 
additionally, calculations are compared to angular dis- 
tributions at hi gher energies and to excitation functions 
from literature [12L HH . Finally, conclusions are drawn 
in Sec. El 



II. EXPERIMENTAL TECHNIQUE 

The experiment was carried out at the cyclotron lab- 
oratory of ATOMKI, Debrecen. A similar experimental 
setup was used also in the previous experiments 0-[H| 
and is described in more detail in [2(J. The following 
paragraphs provide a short description of the experimen- 
tal procedure. 



A. Targets and scattering chamber 

The targets were produced by evaporation of metallic 
yttrium onto thin carbon foils (« 20 /ig/cm 2 ). The target 
thickness was approximately 200 /jg/cm 2 . The targets 
were mounted on a remotely controlled target ladder in 
the center of the scattering chamber. The stability of the 
targets was checked continuously with monitor detectors 
(see below) during the experiment. 

The E lab = 16.21 and 19.47 MeV energy alpha beam 
(beam current: 150 pnA) was lead to the scattering 
chamber through an analyzing magnet. Since the energy 
stability of the beam is important for the experiment, 
the widths of the slits at the entrance and the exit of the 
magnet were small (1 mm). Moreover, the beam energy 
was monitored during the experiment with the monitor 
detectors (see below) . The total uncertainty of the beam 
energy was found to be less than 0.5%. 

An aperture of 2 x 6 mm was mounted on the tar- 
get ladder to check the beam position and size of the 
beamspot before and after every change of the beam en- 
ergy or current. We optimized the beam until not more 
than 1% of the total beam current could be measured on 
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this aperture. As a result of the procedure, the horizontal 
size of the beamspot was below 2 mm during the whole 
experiment which is crucial for the precise determination 
of the scattering angle. 

B. Detectors and data acquisition 

Four ion implanted silicon detectors with active areas 
of 50 mm were used for the measurement of the angu- 
lar distributions. The detectors were mounted in pairs 
separated by 10°. The solid angles covered by the two 
detector pairs were Af2=1.56 x 10~ 4 and Af2=1.81 x 
10~ 4 . The ratios of solid angles of the different detec- 
tors were checked by measurements at overlapping angles 
with good statistics. 

In addition, two detectors were mounted at a larger 
distance on the wall of the scattering chamber at fixed 
angles z9=±15° left and right to the beam axis. These 
detectors were used as monitor detectors during the ex- 
periment to normalize the measured angular distribution 
and to determine the precise position of the beam on the 
target. The solid angle of these detectors was Af2=8.2 x 
10" 6 . 

The signals from all detectors were processed using 
charge-sensitive preamplifiers. The output signals were 
further amplified by main amplifiers and fed into analog- 
to-digital converters. Since the elastic scattering cross 
sections at forward angles differ several orders of mag- 
nitude from the one measured at backward angles, a re- 
liable dead time correction is crucial. The data were 
collected using the WinTMCA system which provides an 
automatic dead time control. This automatically deter- 
mined dead time was verified using a pulser in all spectra. 

The energy of the first of the excited state of the 89 Y 
nucleus is 908.97 keV [2l|. There is a large difference 
between the spin of the ground and the first excited states 
(1/2 - and 9/2+ respectively). Therefore the expected 
inelastic scattering cross section is very low (below 10~ 4 
mbarn, calculated with the DWUCK code [22|]) at the 
measured energies. 

This fact explains why events corresponding to inelas- 
tic alpha scattering on 89 Y are missing from the spectra. 
Typical spectra are shown in Fig. [1] The relevant peaks 
from elastic 89 Y-a scattering are well separated from elas- 
tic and inelastic peaks of target contaminations, and - as 
expected - peaks from inelastic a scattering on 89 Y are 
not visible. 



C. Angular calibration 

Knowledge on the exact angular position of the de- 
tectors is of crucial importance for the precision of a 
scattering experiment since the Rutherford cross section 
depends sensitively on the angle. The uncertainty in the 
angular distribution is dominated by the error of the scat- 
tering angles in the forward region. A tiny uncertainty of 



A-d = 0.3° results in a significant error of approximately 
5% in the Rutherford normalized cross sections at very 
forward angles. 

To determine the scattering angle precisely, we mea- 
sured kinematic coincidences between elastically scat- 
tered alpha particles and the corresponding 12 C recoil 
nuclei using a pure carbon foil target. One detector was 
placed at $ = 70°, and the signals from the elastically 
scattered alpha particles on 12 C were selected as gates 
for the other detector which moved around the expected 
12 C recoil angle $ = 45.85°. Figure [5] shows the relative 
yield of the 12 C recoil nuclei in coincidence with elasti- 
cally scattered alpha particles as a function of the 12 C 
recoil angle. The final angular uncertainty was found to 
be Ai? < 0.15°. 



D. Elastic scattering cross sections and data 
analysis 

Complete angular distributions between 20° and 170° 
were measured at energies of E Q = 16.21 and 19.47 MeV 
in 1° (20° < i? < 100°), 1.5° (100° < d < 140°) and 2° 
(140° < d < 170°) steps. 

The statistical uncertainties varied between 0.1% (for- 
ward angles) and 4% (backward angles). The count rates 
Nft) have been normalized to the yield of the monitor 
detectors N Mon .(tf=15°): 

\dn) [ ' \dn) Mon Nmou. An ' u 

with AJ7 being the solid angles of the detectors. The 
cross section at the position of the monitor detectors 
is taken as pure Rutherford. The relative measurement 
eliminates the typical uncertainties of absolute measure- 
ments, coming mainly from changes in target and from 
the beam current integration. 

The measured angular distributions are shown in 
Fig. [5J The lines are the result of optical model predic- 
tions using global a- nucleus potentials (see Sec.lIVp. The 
measured absolute cross sections cover five orders of mag- 
nitude between the highest (forward angles at E Q =16.21 
MeV) and the lowest cross sections (backward angle at 
E Q =19.47 MeV) with almost the same accuracy (4-5% 
total uncertainty). This error is mainly caused by the 
uncertainty of the determination of the scattering angle 
in the forward region and from the statistical uncertainty 
in the backward region. 

III. LOCAL OPTICAL POTENTIAL AND 
a-CLUSTER STATES 

A. Local folding potential 

The complex optical model potential (OMP) is given 
by: 
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FIG. 1: Typical spectra at E c . m . = 15.51 MeV and 18.63 MeV, # c . m . » 40° and 150°. The peak from elastic 89 Y-a scattering 
is well resolved from both the 12 C-a and 16 0-a elastic scattering. The pulser peak used for the dead time correction is also 
shown. Note the logarithmic scale of all snectra. 
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FIG. 2: Relative yield of 12 C recoil nuclei in coincidence with 
elastically scattered a particles. The shaded area represents 
the angle and uncertainties expected from the reaction kine- 
matics. A Gaussian curve fitted to the experimental data is 
shown to guide the eye. 



U(r) = V c (r) + V(r) + iW(r), 



(2) 



where Vc(r) is the Coulomb potential, V(r), and W(r) 
are the real and the imaginary parts of the nuclear po- 
tential, respectively. The real part of the potential is 
calculated from the folding procedure [H, 0] using a 
density-dependent nucleon-nucleon interaction. The cal- 



culated folding potential is adjusted to the experimental 
scattering data by two parameters 



V(r) =\V F (r/w) 



(3) 



where A ~ 1.1 — 1.4 is the potential strength param- 
eter Q and w « 1.0 ± 0.04 is the width parameter that 
slightly modifies the potential width. (Larger deviations 
of the width parameter w from unity would indicate a 
failure of the folding potential.) The nuclear densities of 
89 Y and a are derived from the measured charge density 
distributions [25| . For details of the folding potential see 
also [IE [13. 

The imaginary part W(r) is taken in the usual Woods- 
Saxon parametrization. For the fits to the experimental 
data we use a sum of volume and surface potential: 



W(r) = W v x f(x v ) + 4 W s x 



df(xs) 
dxs 



(4) 



with the potential depths Wv and Ws of the volume 
and surface parts and 



f{Xi) = 



1 



1 + cxp (Xi 



(•5) 



1 /3 

and Xi = (r — Ri * Aj, ) /aj with the radius parameters 
Ri, the diffuseness parameters a^, and i = S, V. 
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FIG. 3: Rutherford normalized elastic scattering cross sec- 
tions of 89 Y(a, a) 89 Y reaction at E c . m . = 15.51 and 18.63 
MeV versus the angle in center-of-mass frame. The lines are 
the results from the local potential which is fitted to the ex- 
perimental data. The parameters of the fits are listed in Table 
11 



The adjustment of the OMP parameters leads to an 
excellent description of the new experimental data (see 
Fig. [3]). The parameters of the potentials are listed in 
Table H 

It has to be noted that the folding potential in the real 
part may be replaced by a Woods-Saxon potential. A 
similar quality of the fits can be obtained in this case. 
However, the adjustment of three Woods-Saxon parame- 
ters in the real part and six parameters in the imaginary 
part leads to a variety of parameter sets with compara- 
ble x 2 - This problem has been reduced by [3] where a 
microscopic Woods-Saxon potential was derived from a 
folding procedure. This potential is used as the base of 
one of the global potentials studied in Sect. IIVI 

The above calculation with folding potentials has been 
repeated for the 89 Y(a,a) 89 Y scattering data available 
in literature from about 20 MeV up to about 170 MeV. 
The analyzed data have been measured by Brissaud et 
al. HI at 166 MeV (E c . m . = 158.9 MeV), Bingham et 
al. [13 at 65 MeV (E c . m . = 62.2 MeV), Alster et al. 
at 42 MeV (E c . m . = 40.2 MeV), England et al. [3 at 
25 MeV (E c . m , = 23.9 MeV), and Wit et al. [13 at 25, 
23, and 21 MeV (E c . m . = 23.9, 22.0, and 20.1 MeV) 
with an additional excitation function at backward an- 
gles. The data measured by Berinde et al. [17[ at 27.3 
MeV (-Ecm. = 26.1 MeV) were excluded because all fits 
were of poor quality and required width parameters w 
deviating strongly from unity. The new and challenging 
technique for the measurement of angular distributions 
using a position-sensitive detector in the work of Berinde 
et al. [171 ] was unfortunately not further developed; this 
may indicate that there are problems with the data of 



0- 

At 25 MeV two data sets are available. The data by 
England et al. [13] cover the whole angular range from 
forward to backward directions. As expected, at forward 
angles the data agree with the Rutherford cross section, 
and thus the uncertainty of the absolute normalization of 
the data is of the order of a few per cent. The data by Wit 
et al. [l2| focus on the backward region and start around 
"#c.m. ~ 50° where the cross section is about 20 % — 30 % 
of the Rutherford cross section. Although an uncertainty 
of about ±15 % is claimed in [12j, the comparison to the 
data by England et al. shows that the data by Wit 
et al. have to be reduced by a factor of 1.45 to come 
into agreement with the England et al. data. The ori- 
gin of this discrepancy remains unclear. Because typi- 
cal normalization problems (e.g. uncertainties of the tar- 
get thickness) apply probably to all measurements of the 
89 Y(a,a) 89 Y cross section in [l2j, we use the reduction 
factor of 1.45 for all data for 89 Y measured by Wit et 
al. [HI], i-e. the three angular distributions at 21, 23, and 
25 MeV and the excitation function at backward angles. 

The calculated cross sections are compared to the ex- 
perimental data in the broad energy range from about 
20 MeV up to about 160 MeV in Fig. H Similar to the 
result for the new low-energy data shown in Fig. [31 excel- 
lent agreement is obtained for all angular distributions in 
the wide energy range from slightly above the Coulomb 
barrier up to 160 MeV. The parameters of the potential 
are also listed in Table HI 

Unfortunately, the angular distributions at E c , m , = 
40.2 MeV and 62.2 MeV do not cover the backward an- 
gular range which is most sensitive to the imaginary part 
of the potential. On the other hand, the derived poten- 
tials from the analysis of the England et al. data at 25 
MeV and the corrected Wit et al. data at 25 MeV are 
in reasonable agreement; this indicates that the Wit et 
al. data that focus on the backward angular region are 
sufficient to derive the optical potential - provided that 
the absolute normalization is correct. In practice this 
means that data in the forward region are not necessary 
to derive the potential but are highly necessary to define 
the absolute normalization of the data by comparison to 
the Rutherford cross section. 

The various available scattering data for 89 Y enable a 
study of the energy dependence of the potential param- 
eters that are derived from the fits to the angular distri- 
butions in Figs. [3] and |4] For the real part a very smooth 
dependence of the strength parameter A and the width 
parameter w is found, see Fig.[S] As expected, the width 
parameter w remains close to unity at all energies, and 
the strength parameter A is almost energy-independent 
and varies less than 10% between 1.27 and 1.35 leading 
to volume integrals Jr for the real part in agreement with 
the systematics of 0] • The much stronger decrease of the 
real volume integral Jr with energy by about 25 % is a 
consequence of the energy dependence of the interaction 
and the dispersion relation which couples the real and 
imaginary parts of the potential. (Note that - as usual 
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TABLE I: Parameters of the local potentials that were derived from elastic scattering angular distributions in a wide energy 
range (see Figs. O and Q. 



E c .m. 


A 


w 


Jr 


fR,rms 


Ji 


fl ,rms 


Wvo 


Rv 


ay 


4 x Wso 


Rs 


as 


experimental 


(MeV) 






(MeVfm 3 ) 


(fm) 


(MeVfm 3 ) 


(fm) 


(MeV) 


(fm) 


(fm) 


(MeV) 


(fm) 


(fm) 


data from 


15.5 


1.350 


0.986 


342.3 


4.905 


49.9 


5.953 


-11.3 


1.656 


0.486 


-8.3 


1.350 


0.548 


this work 


18.6 


1.353 


0.979 


335.3 


4.868 


46.9 


4.975 


-20.0 


1.331 


0.713 


-25.8 


1.487 0.255 


this work 


20.1 


1.328 


0.985 


334.9 


4.897 


45.0 


5.262 


-22.3 


1.197 


0.982 


-55.7 


1.501 


0.174 


Ref. [12] 


22.0 


1.352 


0.969 


324.7 


4.817 


47.2 


5.856 


-16.9 


1.585 


0.822 


-42.7 


1.460 


0.502 


Ref. [12] 


23.9 


1.315 


0.989 


335.8 


4.917 


46.6 


4.955 


-21.9 


1.656 


0.590 


-71.2 


1.435 


0.605 


Ref. [12] 


23.9 


1.326 


0.984 


334.5 


4.896 


44.7 


4.735 


-23.6 


1.672 


0.610 


-78.5 


1.423 


0.678 


Ref. [13] 


40.2 


1.268 


1.002 


328.0 


4.988 


57.5 


5.672 


-18.8 


1.492 


0.799 


-26.3 


1.475 


0.429 


Ref. [16] 


62.2 


1.296 


1.007 


328.3 


5.016 


65.6 


5.835 


-19.5 


1.634 


0.544 


-56.7 


1.388 


0.365 


Ref. [15] 


158.9 


1.317 0.985 


260.0 


4.930 


82.3 


5.785 


-22.1 


1.457 


0.757 








Ref. [14] 



- the negative sign of the volume integrals Jr and Jj is 
neglected in the discussion.) 

The volume integral Jj of the imaginary part increases 
with energy, but it remains below Jj = 100 MeVfm 3 , 
in agreement with the systematics shown in Q. How- 
ever, the slope of J/ vs. E c m is not well-defined from 
the available data, in particular not in the low-energy re- 
gion that is most relevant for nuclear astrophysics. All 
data between 15 MeV < £? c .m. < 25 MeV lead to volume 
integrals of about Jj rj 50 MeVfm 3 . However, a closer 
look at the imaginary potentials shows that the shape 
of the potentials changes significantly between the differ- 
ent energies (see Fig. [SJ lower part) which complicates 
the analysis of the Ji(E) dependence. Contrary to the 
imaginary potentials, the shape of the real potentials is 
well-defined from the folding procedure (see Fig. [6] upper- 
part ) . 

It has to be noted that there are discrete ambiguities 
for the real part of the OMP. This has been illustrated 
e.g. in Fig. 5 of [ll| where 11 so-called "families" of poten- 
tials were identified which resulted in a similar descrip- 
tion of the low-energy 144 Sm(a,a) 144 Sm elastic scatter- 
ing data. However, these discrete ambiguities are sig- 
nificantly reduced in the present analysis because of the 
data at higher energies. It is not possible to describe the 
data at 62.2 MeV and 158.9 MeV using potentials from 
another family, i.e. using a potential strength which is 
increased or decreased by about 30%. 

Together with the determination of the bound state po- 
tential in Sect. lHTBl the family with J R rj 320 - 350 MeV 
fm 3 has been selected in the analysis of the low-energy 
scattering data. For completeness it has to be pointed 
out that the above volume integral of Jr « 320 — 350 
MeVfm 3 at low energies is only valid for folding poten- 
tials. Slightly different numbers will be obtained if other 
parameterizations of the potential, e.g. Woods-Saxon po- 
tentials, are used (see e.g. Q). 

In addition to the angular distributions, an excitation 
function at the ver y b ackward angle of i?i a b = 176° has 
been measured in [12| in the energy range of 18 MeV 
< ^a.iab < 26 MeV, i.e. covering the energ y ra nge of the 
angular distributions in this work and in [121 Il3| . This 



excitation function has been calculated using the poten- 
tial from the fit to the 20. 1 MeV data from [if] , see Fig. [7] 
Additionally, an averaged potential has been calculated 
from all fits below 25 MeV. Both the 20.1 MeV potential 
and the average potential are able to reproduce the shape 
of the excitation function including the deep minimum 
around 23 MeV; however, the agreement between theory 
and experiment is slightly worse for the excitation func- 
tion at very backward angles (compared to the excellent 
reproduction of the angular distributions). Because of 
the relatively weak energy dependence of the volume in- 
tegrals Jr and Ji below 25 MeV, this averaged potential 
will also be used for comparison with global potentials in 
Sect. [TV] 

B. a-cluster states in 93 Nb = 89 Y ® a 

a clustering in intermediate and heavy nuclei is a well- 
established phenomenon [JJ, |28h30| | . A link between the 
cluster model and the collective model has been discussed 
in [3l|, and it has been pointed out that a cluster model 
interpretation of collective properties - as recently re- 
viewed in [32] - emerges. 

The properties of a clustering above shell closures in 
the even-even nuclei 44 Ti = 40 Ca <g> a, 94 Mo = 90 Zr 
<X> a, and 212 Po = 208 Po (g> a have been studied exten- 
sively @, [28T - [3fl I33T438T ] . However, only few studies are 
available for even-odd nuclei below the above mentioned 
shell closures f39rl4i1]. Following the ideas in |H,|4(J for 
43 Sc = 39 K (g) a, we analyze here a-cluster properties of 
93 Nb = 89 Y <g, a . The ground state of 89 Y (Z = 39, 
N = 50) has J T = l/2~ which comes from a one-proton 
hole in the p/1/2 shell of the neighboring 90 Zr (Z = 40, 
N = 50). Counterparts of the well-studied rotational 
bands in 94 Mo = 90 Zr ® a should be found in 93 Nb = 
89 Y ® a with 89 Y = 90 Zr ® p~L. This study extends a re- 
cent review on a-cluster states in N = 50 <g> a even-even 
nuclei [42j . All excitation energies and decay properties 
of excited states in 93 Nb have been taken from [431 ] . 

The full formalism of the applied model can be found in 
e.g. @,[38|. Here we briefly repeat the important features 
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FIG. 5: Potential parameters A and w and the integral po- 
tential strengths Jr and Ji for the real and imaginary part 
of the potentials that are derived from the fits to the angular 
distributions in Figs. [3] and [4] (full circles). Additionally, the 
result from the bound state adjustment in Sect. MI Bl is shown 
with triangles at E = 0. Further discussion see text. 



FIG. 4: Rutherford normalized elastic scattering cross sec- 
tions of 89 Y(a, a) 89 Y reaction at energies between 20 and 160 
MeV. The experimental data are taken from literature 
[l6| . The parameters of the fits are listed in Table [T] 

of the model. The a-cluster wave function can be directly 
calculated from the Schroedinger equation and the a- 
nucleus potential which is taken from the double-folding 
procedure (see also above). The so-called Wildermuth 
condition ensures that the Pauli principle is taken into 
account: 

4 4 

Q = 2N + L = J2{2n i + l i )=J2q i (6) 

i=l i=l 

where Q is the number of oscillator quanta, N is the 
number of nodes and L the relative angular momentum 
of the a-cluster wave function, and qi = 2ni + U are the 
corresponding quantum numbers of the nucleons in the 
a cluster. We take q = 4 and thus Q — 16 for the 
band in 93 Nb = 89 Y(l/2-) ® a. A 1/2+ band with Q = 
17 should also exist in 93 Nb; however, no firm assignment 
could be found up to now [43| ■ Also a higher nodal band 
with Q — 18 has not yet been found experimentally |43j . 
Typical properties of a-cluster states are large reduced 
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FIG. 6: (Color online) Shape of the potentials for the fits in 
Figs.[3]and[31 The upper part shows the very regular behavior 
of the real part V(r). The lower part shows the variations 
of the shape of the imaginary potential W(r) at the various 
energies. Further discussion see text. 
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widths Q\ or spectroscopic factors. Unfortunately, no a 
transfer data like e.g. ( 6 Li,d) on 89 Y can be found in 
43] ; such data have been essential for the assignment of 
a-cluster properties in many cases |44j |. 

A prerequisite for a successful description of a-cluster 
states is an a-nucleus potential which is also able to de- 
scribe the elastic a scatte ring cross section and/or a- 
decay properties (e.g. [H, l38l. |39|). In the present in- 
vestigation of 93 Nb we use a double-folding potential for 
the real part that is close to the result of the previous 
Sect. IIII A[ the imaginary part of the potential vanishes 
at the very low energies studied here. 

The properties of the a-cluster states in 93 Nb with 
a 89 Yb(l/2~) ® a structure have been analyzed in the 
following way. In a first step the strength parameter A 
of the double-folding potential is adjusted to reproduce 
the binding energy of the first 1/2" state in 93 Nb; it is 
located at E x = 30.8 keV or E — Q a -E x = -1900. 7keV 
with the a binding energy Q a = — 1931.5keV. The result 
A = 1.1951 and the corresponding volume integral Jr = 
326.1 MeV fm 3 are very close to the numbers obtained 
from the analysis of the scattering data. For simplicity, 
the width parameter w has been fixed at w = 1 for the 
bound state calculations. 

For a perfect rotator the calculations with a realistic 
potential should be able to predict the energies of all 
members of this rotational band. However, it has been 
shown that the strength of folding potentials has to be 
slightly reduced with increasing angular momentum (or 
increasing excitation energy). It was found that an excel- 
lent description of the energies within rotational bands is 
obtained using a potential strength with a weak depen- 
dence on the angular momentum L: 

X(L) = X(L = 0) - c x L (7) 

with small values for the constant c » (3 — 5) x 10~ 3 
[H,[39[ . Thus, in a second step the potential strength was 
adjusted to the centroid of the 3/2~ and 5/2" states with 
L = 2; we find A = 1.1856 and derive c = 4.75 x 10~ 3 in 
excellent agreement with several other nuclei Then 
we use a weak spin-orbit potential Vls ~l/rx dVp/dr; 



TABLE II: Properties of the Q = 16, K 71 = 1/2" rotational 
band in 93 Nb. For L > 6 only the centroid of the two states 
with J = L ± 1/2 is given. Additionally the excitation energy 
of the 1/2+ , 3/2+ band head of the Q = 17, K* = 1/2+ band 
is predicted. 



l ,r 


Q E 


E x ( 93 Nb) 


A 


— — 


- (keV) 


(keV) 


— 


1/2 - 


16 -1901 


31" 


1.1951" 


2 3/2- 


16 -1244 


687" 


1.1856" 


2 5/2 _ 


16 -1121 


810" 


1.1856" 


4 7/2 _ 


16 -457 


1474 6 


1.1761 6 


4 9/2- 


16 -238 


1693 6 


1.1761 6 


6 11/2-, 13/2- 


16 +613 


2545* 


1.1666 6 


8 15/2 - , 17/2- 


16 +1681 


3613 6 


1.1571 6 


10 19/2-, 21/2" 


16 +2878 


4809 6 


1.1476 6 


12 23/2-, 25/2" 


16 +4221 


6152 6 


1.1381 6 


14 27/2", 29/2" 


16 +5749 


7680 6 


1.1286 6 


16 31/2", 33/2" 


16 +7514 


9445 6 


1.1191 6 


1 1/2+ , 3/2+ 


17 +4765 


6696 b 


1.1904 6 


l/2~ 


18 +9432 


11363 6 


1.1951 6 



a E from experimental data; A adjusted to fit the energy E. 
b E predicted using A(L) = X(L = 0) - c X L, Eq. (0. 

the strength of this spin-orbit potential is adjusted to 
reproduce the splitting of the 3/2" and 5/2" states with 
L = 2. Now all parameters of the potential are fixed, and 
it is possible to calculate the energies of the members of 
the Q = 16 rotational band in 93 Nb. The results are 
listed in Table HU 

There are a number of candidates for the L — 4 states 
with ,r = 7/2- and 9/2" at E x = 1364 keV, 1500 keV, 
1603 keV, and 1916 keV; however, no firm assignment is 
possible at the present stage. Therefore, for the higher- 
lying members of the rotational band with L > 6 only 
the centroid energy of the states with J = L ± 1/2 are 
listed in Table [TT1 In addition, the band heads of the L = 
odd, Q = 17 and L — even, Q — 18 bands are predicted. 

The uncertainties of the predictions in Table [TT| can be 
estimated from the uncertainty of the potential strength 
parameter A(L) and thus from the uncertainty of the con- 
stant c of about 20 % [39( . The uncertainty of the poten- 
tial strength A(L) in Eq. ([7]) increases with L from about 
2 x IO" 3 at L = 4 to about 0.01 around L = 10. The 
resulting uncertainty is of the order of about 150 keV for 
the L = 4 states and up to several hundred keV for states 
with higher angular momentum L. Together with an ad- 
ditional uncertainty from the weak spin-orbit potential 
of the order of 100 keV the total uncertainty for the pre- 
dicted excitation energies is about 200 keV for the L = 4 
states and larger for higher L. Slightly varying potentials 
have been used for the description of bands with different 
Q. This leads to an increased uncertainty of the order 
of 1 MeV for the band heads of the Q = 17 and Q = 18 
rotational bands in 93 Nb. 

Electromagnetic decay properties may provide addi- 
tional hints for the a-cluster structure of the Q = 16 ro- 
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tational band in 93 Nb. It is straightforward to calculate 
the wave functions of these states and the reduced tran- 
sition strengths B(E2, J? — > J J) for transitions from an 
initial state with Jf to a final state with JJ [46| . The re- 
sults are listed in Table IIII1 Indeed, enhanced transition 
strengths of the order of lOW.u. are found for several 
transitions. 

Experimental data for transition strengths are rare and 
have huge uncertainties. Within these uncertainties, the 
transition strengths for the transitions from the L = 2 
states with Jf = 3/2~ and 5/2~ to the final state with 
L = 0, JJ = l/2~ are calculated correctly. 

For transitions from L = 4 states to L = 2 states pre- 
dictions are given in Tabic IIIII It is interesting to note 
that the absolute transition strengths Tt-^f depend sen- 
sitively on the predicted energies of these states because 
of the E^ dependence of E2 transitions. However, the 
reduced transition strengths are almost independent of 
the predicted energy. From the above mentioned can- 
didates for the L = 4 states, only the transition from 
the J t = 7/2 state at E x = 1500 keV to the J] = 5/2" 
state at E x = 810 keV has been detected experimentally 
The preferred electromagnetic decay mode of most of the 
L = 4 candidates proceeds via El transitions to low- 
lying JJ = 7/2+ and 9/2+ states which can be found 
at E x = OkeV, 744 keV, and 1083 keV. Thus, the predic- 
tions of B(E2) values in Table Hill are not able to firmly 
assign the L = 4 members of the Q = 16 rotational 
band because of missing experimental data for weak de- 
cay branches. It has to be noted that a 1 MeV transition 
in 93 Nb with a El strength of 10 -3 W.u. corresponds to 
r 7 (_El) « 1.4meV whereas a similar 1 MeV E2 transi- 
tion with 10 W.u. corresponds to a much smaller radia- 
tion width of only Y 1 (E2) = 0.2 mcV. 

In a recent paper, the state at E x = 1500 keV has 
been assigned J* = 9/2" [H instead of J* = 7/2^ 
[H E3, and an E2 transition strength of 26.4^6.2 W.u. 
has been measured for the transition to the 5/2~ state at 
E x = 810 keV [48]. A tiny readjustment of the potential 
strength parameter A of about 0.2 % from the extrapo- 
lated value of 1.1761 to 1.1788 is required to shift the 
9/2~ state from the predicted value E x = 1693 keV to 
1500 keV. Simultaneously, the 7/2~ state is shifted from 
the prediction E x = 1474 keV to 1282 keV which is very 
close to the J w = 5/2~, 7/2" state at E x = 1364keV. 
As pointed out above, the calculated transition strength 
of 12 W.u. for the 9/2~ — > 5/2~ E2 transition is practi- 
cally independent of the precise excitation energies. The 
enhanced experimental transition strength confirms the 
,r = 9/2- assignment of [H for the E x = 1500 keV 
state and the interpretation as a member of the L = 4 
doublet with 89 Y <g> a structure. 



IV. GLOBAL OPTICAL MODEL PREDICTIONS 

After the successful study of 89 Y(a,a) 89 Y elastic scat- 
tering and 93 Nb = 89 Y ® a bound state properties with 



a locally optimized potential we will now test the pre- 
dictions of global a-nucleus optical potentials for the 
89 Y(a,a) 89 Y elastic scattering cross section. Obviously, 
the best description of elastic scattering data is obtained 
from a locally adjusted potential. However, such a lo- 
cally optimized potential requires experimental data for 
adjustment which are often not available. In particular, 
for unstable nuclei it is not possible to measure elastic 
scattering data with the required precision. Therefore 
global a-nucleus potentials are required for the predic- 
tion of scattering and reaction data. In the following we 
will compare the predictions of published global poten- 
tials to our experimental results. 

The prediction of cross sections for unstable nuclei re- 
quires the extrapolation of global potentials into regions 
of the nuclear chart where no scattering data are avail- 
able to test the extrapolations. It has been shown re- 
cently [l(| that the variation of the potential along an 
isotopic chain (in that case Z = 50) can be sensitively 
studied by a comparison of the scattering cross sections 
of a neutron-rich and a neutron-deficient isotope. Here 
we extend this idea and study the variation of the scat- 
tering cross sections along the N = 50 isotonic chain by 
a comparison of scattering data for 89 Y and 92 Mo. 

Several different parameterizations for the optical po- 
tential exist. For reactions involving alpha particles the 
optical potential is described mainly by complex Woods- 
Saxon potentials. Some authors investigated the use of 
hig her order terms of the Woods-Saxon function, too 
jig . [50j . Mo del- independent parameterizations have also 
been studied, either with spline functions Ht|[5lJ or series 
of Fourier-Bessel functions added to the Woods-Saxon 
parameterizations [521 [53[ or with a sum of Gaussians 
[50| . A folding model has been introduced by Kobos et 
al. [54j . In the framework of this model an energy- and 
density-dependent effective interaction — approximating 
the nonlocal exchange component by an empirical param- 
eterization — was used to describe alpha-particle elastic 
scattering at energies from 25 to 120 MeV. This model 
was also extended for inelastic scattering. However, in all 
of these investigations the parameterization of the imag- 
inary part of the optical potential varied from nucleus 
to nucleus, in order to obtain a good description of the 
scattering cross sections for all the nuclei studied. 

In the framework of the p process network calcula- 
tions a large number of reactions involving alpha particles 
(alpha-induced reactions and alpha particle emission) has 
to be taken into account. As the p process path is located 
in a region of unstable nuclei on the neutron-deficient 
side of the chart of nuclides, experimental data are prac- 
tically not available to adjust potential parameters of the 
a-nucleus potential. Therefore, a global a-nucleus opti- 
cal potential is required for the theoretical prediction of 
reaction cross sections involving a particles within the 
statistical Hauser-Feshbach model. 

Several different global and regional parameterizations 
have been developed in recent years to describe the inter- 
action between nuclei and alpha particles. In the follow- 
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TABLE III: Calculated reduced transition strengths B(E2) in the Q = 16, K n — 1/2 rotational band in 93 Nb compared to 
experimental values [43| . 
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B(E2, Ji — » J/) ca i c 
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(keV) 
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3/2" 
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1/2" 


31 


8.6 


1J --10 


2 
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810 





1/2" 


31 


8.6 


< 79 


4 


7/2- 


1474 


2 


3/2- 


687 


10.7 




4 


7/2- 


1474 


2 


5/2- 


810 


1.2 




4 


9/2- 


1693 


2 


5/2- 


810 


12.0 





ing an overview of the parameterizations is given which 
are studied in this work. 

(i) The regional a-nucleus potential of Avrigeanu 
ct al. [H corresponds to average mass-, charge-, and 
energy-dependent Woods-Saxon parameters based on lo- 
cal potentials obtained by analysis of 108 experimental 
angular distributions of a-particle elastic-scattering on 
target nuclei from 50 Ti to 124 Sn and a-particle energies 
from 8.1 to 49MeV. The local Woods-Saxon parameter 
sets provided by the analysis of the new experimental 
data (Table 2 of Ref. 0) have led to similar results as 
those shown in Fig. [3] (see also Fig. 7 of Ref. Q). 

(m) The recently published global potential by Ku- 
mar et al. Q claims to describe data in a wide mass 
(12 < A < 209) and energy region (from the Coulomb 
barrier up to 140 MeV); however, problems below 30 MeV 
had to be compensated by an enhanced imaginary volume 
integral, and even with this enhancement the description 
of 90 Zr(a,a) 90 Zr scattering at 15 MeV is reasonable but 
not excellent (see Fig. 2 in [5]). A similar reproduction 
of the 89 Y(a,a) 89 Y scattering data below 20 MeV is ex- 
pected. 

(Hi) The widely used potential by McFadden and 
Satchler [59| is a very simple 4-parameter Woods-Saxon 
potential with mass- and energy-independent parame- 
ters. Despite its simplicity it provides an excellent de- 
scription of a-scattering data and cross sections of a- 
induced reactions; e.g., this potential is used as default 
for the Hauser-Feshbach calculations of astro phy sical re- 
action rates by Rauscher and Thielemann [3, |55|. 

Unfortunately, the latest version (5(| H?} of the poten- 
tials by Demetriou et al. Q are only published in confer- 
ence proceedings and cannot be used without the authors 
of [56|,[57|; we do not intend to show results from the early 
and perhaps out-dated potentials in Q. 



A. Angular distributions: comparison with 
theoretical models 

The regional optical potential (ROP) of Avrigeanu et 
al. was obtained by a semi-microscopic analysis, us- 
ing the double folding model [58[ , based on alpha-particle 
elastic scattering on A w 100 nuclei at energies below 
32 MeV. The energy-dependent phenomenological imagi- 
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FIG. 8: (color online). Rutherford normalized elastic scatter- 
ing cross sections of 89 Y(a, a) 89 Y reaction at E c .,„. = 15.51 
and 18.63 MeV versus the angle in center-of-mass frame. 
The lines correspond to the predictions using different global 
optical potential parameter sets. For more information see 
SeclWl 



nary part of this semi-microscopic optical potential takes 
into account also a dispersive correction to the micro- 
scopic real potential. A small revision of this ROP and 
especially the use of local parameter sets were able to 
describe the variation of the elastic scattering cross sec- 
tions along the Sn isotopic chain [18j. A further step 
to include all available a-induced reaction cross sections 
below the Coulomb barrier has recently been carried out 
First, the ROP based entirely on a particle elastic 
scattering Q was extended to A ~ 50-120 nuclei and en- 
ergies from ~ 13 to 50 MeV. Secondly, an assessment of 
available (a, 7), (a,n) and (a,p) reaction cross sections 
on target nuclei ranging from 45 Sc to 118 Sn at incident 
energies below 12 MeV was carried out. In this work the 
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TABLE IV: Xred °f predictions using diflerent global param- 
eterizations compared with the angular distributions derived 
in the present work. No parameters have been adjusted to 
the new experimental data. 



Global Parameterizations 


89 Y(a, 


a) 89 Y 


Ref. 




15.51 MeV 


18.63 MeV 




Avrigeanu 


6.5 


5.5 


f4] 


Kumar 


181 


287 


[5] 


McFaddden 


35.3 


40.9 


[59] 



10% 




89w, >89w 

Y(a,a) Y 






E c m = 22.0 MeV 


10' : 






X 






t 












b 






10 S 


■ Wit etal., corrected 
- Avrigeanu etal. 
Kumar etal. 




10"S 


McFadden era/. 




10°- 




Y(a,a) Y 




E c m = 23.9 MeV 


10" 1 - 






X 






h- 






Z> 

tr 


\ J' ji- 




e io 3 - 













10" 3 - 


■ England ef a/. 

Avrigeanu era/. 

Kumar et al. 




10" 4 - 


McFadden etal. 


V V ¥ 



20 40 60 80 100 120 140 160 180 



f> [deal 

cm. L OJ 

FIG. 9: (color online). Rutherford normalized elastic scatter- 
ing cross sections of 89 Y(a, a) 89 Y reaction at E c . m . = 20.1, 
22.4 and 23.9 MeV versus the angle in center-of-mass frame 
taken from Wit et al. and England et al. p3, . The data 
by Wit et al. [12] have been corrected by a factor of 1.45 to 
come into agreement with the England et al. data ,13] at 23.9 
MeV. The lines correspond to the predictions using different 
global optical potential parameter sets. For more information 
see Sec. HV1 



most recent potential of Avrigeanu et al. [4(, is used to 
calculate the scattering cross sections. 

The potential from Kumar et al. Q was proposed to 
describe alpha-induced reactions on 12 < A < 209 target 
nuclei at E < 140 MeV. In that work, the systematics 
of volume integrals has been used to determine the real 
and the imaginary parts of the potential. The real po- 
tential volume integrals have been taken from the work 
of Atzrott et al. [Tj . The best fit volume integrals from 
the phenomenological analysis, consistent with the above 



set, had been employed in their analysis. A similar ap- 
proach was followed for the imaginary part starting with 
the volume integral systematics of the imaginary part. 
Moreover at energies below 30 MeV a dispersive correc- 
tion is taken into account. Calculations using this poten- 
tial are able to reproduce the measured alpha scattering 
angular distributions at higher energies Q . Here we com- 
pare the predicted angular distributions from the Kumar 
et al. [E| potential to our new experimental data at lower 
energies close to the Coulomb barrier. 

For completeness, we take into account also the po- 
tential of McFadden et al. [59j]. Numerous alpha elastic 
scattering experiments were done on nuclei between O 
and U at alpha energies of 24.7 MeV in the 60's. Fits 
were obtained using a four-parameter Woods-Saxon po- 
tential. 

In Figure [5] the measured angular distributions of the 
89 Y(cv, a) 89 Y elastic alpha scattering at E c . m .= 15.51 and 
E c . m .=18.63 MeV are shown. The different lines corre- 
spond to the predictions using the above discussed global 
and regional optical potential parameterizations without 
any further adjustment of parameters. The overall agree- 
ment between the calculations performed with using the 
potentials of Avrigeanu and McFadden and the exper- 
imental data is good. Only the calculation performed 
with the potential of Kumar et al. gives a significantly 
worse description of the experimental data. For a strict 
comparison between the potentials of Avrigeanu, and Mc- 
Fadden a x 2 analysis has been done. The resulting x 2 
parameters can be found in Table I. As can be seen, the 
quality of the different parameterizations is similar al- 
though a value of x 2 e d ~ 1 cannot be reached by any of 
the global potentials. 

Our analysis is extended up to 23.9 MeV using the 
data of Wit et al. [n| and England et al. [H|. The 
89 Y(a, a) 89 Y elastic scattering was studied by Wit et al. 
at E c . m . fa 20.1, 22.4, and 23.9 MeV and by England 
et al. at E c . m . pa 23.9 MeV. Note that we have applied 
the same reduction of a factor of 1.45 to all experimental 
data on 89 Y in the Wit et al. paper [3 (see Sect. HUT) . 

In general, both the Avrigeanu and McFadden parame- 
terizations can describe roughly both the magnitude and 
the oscillation pattern of the angular distributions mea- 
sured at energies between 20.1 and 23.9 MeV. The calcu- 
lations performed with the Kumar potential overestimate 
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FIG. 10: (color online). Excitation functions for elastic a 
scattering on 89 Y target nuclei at a scattering angle of 176° 
taken from the literature [13 ] and corrected by a factor of 1.45 
(see text in Sect. [TV^|) . The different curves correspond to 
the predictions using different global and regional potentials. 
For more information see Sec. IIVI 



the strength of the oscillation, predicting deeper minima 
at backward angles. In the case of calculations performed 
with the other two global parameterizations slight dif- 
ferences at backward angles between the measured and 
calculated data can be seen in Fig. |H1 
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FIG. 11: (Color online) Comparison between the averaged 
local potential (see Sect. nTTX|) and the global parameteriza- 
tions at 18.63 MeV (see Sect HV) . 



ward angles (in particular in the pronounced minimum 
around 23 MeV in the excitation function, see Fi g. [3 
may be affected by small compound contributions [60j . 
Thus, minor deviations between the experimental and 
calculated cross sections are acceptable. 



Excitation functions 



A perfect global optical potential should be able to 
predict angular distributions including the diffraction- 
like patterns which show up at backward angles with 
increasing energy. It has to be pointed out that the 
backward angular region has a higher sensitivity to the 
optical potential than the forward angular region where 
the Coulomb interaction is dominating. For a particular 
study of this backward angular region, an excitation func- 
tion has been measured by Wit et al. [HJ at $ = 176°. 

The measured excitation function of alpha particle 
elastic scattering on 89 Y at a scattering angle 176° [l2j 
and the corresponding optical model calculations using 
the potential from Avrigeanu et al. Kumar et al. [5| 
and McFadden and Satchler [59[ are shown in Fig. [TO] 
Although there is reasonable agreement with the magni- 
tude of the cross section in all calculations, the global po- 
tentials are not able to reproduce the energy dependence 
of the cross section. The sharp minimum of the experi- 
mental excitation function around 23 MeV is rather well 
described by the potential of Avrigeanu et al. Q, sug- 
gested by the global potential of McFadden and Satchler, 
and not predicted by the Kumar et al. global potential. 
The appropriate description of particular features of the 
experimental elastic scattering of alpha particles on 89 Y 
proves conclusively the suitable account of the nuclear 
absorption merely by the potential of Avrigeanu et al. 
However, the small scattering cross section at very back- 



C. Comparison of global and local potentials 

For further improvement of the global a- nucleus poten- 
tials a deeper understanding of the differences between 
the potentials and the resulting cross sections is required. 
In the following paragraphs we will compare the shapes 
of the different potentials and the elastic phase shifts Sl 
and scattered wave amplitudes t)l that define the elastic 
scattering cross section. Note that the scattering cross 
section is related to the sum of the absolute squares over 
all contributing partial waves; the study of the underly- 
ing Sl and rj^ will allow a better understanding of the 
global potentials. However, it is beyond the scope of the 
present paper to derive a new and improved global a- 
nucleus potential. 

The shapes of the real and imaginary part of the local 
potential have already been shown in Fig. [51 Because 
of the minor energy dependence of the volume integrals 
Jr and J i below 30 MeV (see Fig. [5]), an averaged local 
potential is used for comparison to the global potentials. 
The global potentials are shown in Fig. QTJ Huge differ- 
ences can be seen for small radii, i.e. the nuclear interior, 
whereas at larger radii (r > 8 fm) all potentials are close 
to each other. 

In a semi-classical picture larger radii in the potential 
correspond to larger impact parameters that are related 
to the angular momentum L — fx p. E.g., for the energy 
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FIG. 12: (Color online) Comparison between the local po- 
tential (see Sect. IHlX| and the global parameterizations (see 
Sect lIV|l for E = 18.63 MeV: elastic scattering phase shift Sl 
(upper) and amplitude of the scattered wave t)l (lower) in 
dependence of the angular momentum L. The data points 
are connected by dashed lines to guide the eye. Note that the 
scale in L with < L < 22 corresponds roughly the scale in 
r of Fig. ITT1 with < r < 12 fm according to Eq. ©. 



E = 18.63 MeV this leads to the approximate relation 
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FIG. 13: (color online). Rutherford normalized elastic scat- 
tering cross sections of 92 Mo(a, a) 92 Mo reaction at E c . m . = 
15.69 and 18.63 MeV versus the angle in center-of-mass frame 
taken from Fiilop et al Q. The lines correspond to the 
predictions using different global optical potential parameter 
sets. For more information see Sec. IIVI 



\L\/h » 1.85 x r/fm (8) 

As expected, the scattering phase shifts Sl and scattered 
wave amplitudes rjL arc in agreement for L <J 15 (see 
Fig. [T2")) . However, the nuclear potential is relatively 
weak at such large radii, and thus the phase shifts Sl are 
small and t\l ~ 1 (no absorption); the scattering cross 
section is dominated by the Coulomb interaction. 

The situation changes for smaller angular momenta. 
Around L « 10 (or semi-classically, r « 5.4 fm) the de- 
creasing rjL indicate increasing absorption. Simultane- 
ously, the Sl start to deviate from zero. For small L < 5 
(corresponding r < 2.7 fm) the tjl values are coming close 
to zero, i.e. full absorption of the respective partial wave. 
For these partial waves the calculated phase shift Sl does 
not have strong impact on the scattering cross section. 
In other words, the elastic scattering cross section is not 
very sensitive to the potential at small radii. However, 
one should keep in mind that this simple semi-classical 
interpretation is not strictly valid; the determination of 
Sl and t]l from the solution of the Schrodinger equation 
depends on the underlying real and imaginary potentials 
V(r) and W(r) for the whole range in r where V(r) or 
W(r) deviate from zero. 

From Fig. [T^] it can be seen that the potential by 
Avrigeanu et al. provides Sl and tjl that are close to 
the result of the local analysis (which is taken as a ref- 
erence here) for all L. In particular, the most relevant 
partial waves around L « 10 are nicely reproduced. The 



McFadden/Satchlcr potential leads to slightly but sys- 
tematically larger values of rjL for the most relevant L 
around L w 10. This weaker absorption of the McFad- 
den/Satchler potential can also be seen in Figs. [8] and [9] 
where the calculated cross section in the backward region 
overestimates the experimental data. The phase shifts for 
small L are best reproduced by McFadden/Satchler; how- 
ever, as pointed out above, the scattering cross section is 
not very sensitive to this region. (This may be different 
for the calculation of a-induced reaction cross sections 
where the McFadden/Satchler potential has been used 
very successfully!) The predictions of Sl and tjl from 
the Kumar et al. potential deviate significantly from the 
results of the local potential, and consequently the calcu- 
lated scattering cross sections are not in good agreement 
with the experimental data. 



D. Variation of the scattering cross section along N 
= 50 isotonic chain 

Since modeling explosive nucleosynthesis scenarios re- 
quires reaction rates on large number of reactions involv- 
ing thousands of nuclei, the ct-nucleus potential has to be 
known in a wide region. The reliability of the extrapo- 
lation to unstable nuclei can be tested by measuring the 
elastic scattering cross sections on several nuclei along 
isotopic and isotonic chains. The ratio of Rutherford nor- 
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malized cross sections along isotopic or isotonic chains is 
a very sensitive observable for the quality of a-nucleus 
potentials that was not taken into account in most pre- 
vious studies. This interesting feature was not observed 
in earlier work because the backward angular range was 
usually not measured with sufficient accuracy. 

The 92 Mo(a, a) 92 Mo reaction has been investigated by 
Fiilop et al. at E C . OT =13.20, 15.69 and 18.62 MeV Q. 
For completeness, in Fig. Q2] the angular distributions of 
clastically scattered alpha particles on 92 Mo measured at 
E c . m . = 15.69, 18.62 MeV are shown. The different lines 
correspond to theoretical cross sections calculated from 
the above discussed global alpha nucleus potentials, again 
without any further adjustment of parameters. As can be 
seen, the situation is similar to the one found in the case 
of 89 Y(a, a) 89 Y. Namely, the Avrigeanu and McFadden 
global parameterizations can describe the measured an- 
gular distributions with similar quality and the potential 
of Kumar fails to reproduce the experimental data. 

In order to investigate the behavior of the optical po- 
tential parameters along the N = 50 isotonic chain we 
derived the ratio of the elastic scattering cross sections 
of 92 Mo and 89 Y which arc both neutron-magic nuclei. 
It is found that the normalized elastic alpha scattering 
cross sections of 89 Y and 92 Mo differ by roughly 50-70% 
at backward angles, and the ratio shows a pronounced 
oscillation pattern. The large number of experimental 
points and the low uncertainties on both data sets pro- 
vide a unique probe to understand the evolution of the 
a-nucleus potential along the N — 50 isotonic chain. 

In Figure [T3] the experimental ratio of the Rutherford 
normalized elastic scattering cross sections is compared 
to the corresponding results of the recent potential of 
Avrigeanu et al. [H, the most recent global OMP from 
Kumar et al. [f| and the well-known potential of Mc- 
Fadden and Satchler [59(. It can be clearly seen that no 
global parameterization can describe correctly the am- 
plitude and the phase of the oscillation pattern of the 
experimental data at backward angles. This new observ- 
able may provide constraints for the further improvement 
of global a-nucleus potentials. 

It may be added that, while the ROP [4] repro- 
duces only the phase of oscillations rather well, the local 
Woods-Saxon parameter sets corresponding to these tar- 
get nuclei and incident energies (Table 2 of Ref. 
provide a good description of this ratio at 15 MeV and 
a reasonable one at 19 MeV. Therefore, taking into ac- 
count the ROP averaging nature, one may conclude that 
the difference between the experimental ratio values and 
those provided by the ROP describes the variance for 
particular nuclei with respect to the average behavior. 

Recently, a similar study has been performed by 
Galaviz et al. [l(| where the variation of the elastic scat- 
tering cross sections along the tin isotopic chain had 
been studied. Complete angular distributions of the 
112 - 124 Sn(a, a) 112 ' 124 Sn reactions at 18.8 MeV were mea- 
sured. It was found that the elastic alpha scattering cross 
sections of the 112 Sn and 124 Sn differ by roughly 30-40% 
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FIG. 14: (color online). Experimental ratio of the scattering 
cross sections (a/anuTH) 92 Mo / (o-/o-ruth) 89 Y at Ec.m. ~ 
15.51 and 18.63 MeV versus the angle in center-of-mass frame. 
The cross sections of the 92 Mo(a, a) 92 Mo are taken from []|. 
The lines correspond to the predictions using different global 
optical potential parameter sets. For more information see 
Sec. EVJ 

at backward angles, and the ratio of the measured cross 
sections shows a similar oscillation feature. It became 
evident that the global alpha-nucleus potentials failed to 
reproduce either the strength or the oscillation pattern 
for backward angles (Tol j. This behavior is very similar 
to the ratio of the Rutherford normalized cross sections 
of the 92 Mo(a, a) 92 Mo and 89 Y(a, a) 89 Y derived in the 
present work. This fact clearly indicates that the avail- 
able theoretical alpha nucleus optical potential parame- 
terizations have to be improved to be able to describe 
the variation of the angular distributions simultaneously 
along isotopic and isotonic chains. This is particularly 
important for the extrapolation into regions of the chat 
of nuclides where no scattering data exist. 



V. SUMMARY 

In the present work angular distributions of clastically 
scattered alpha particles on 89 Y have been measured at 
E c . m . = 15.51 MeV and 18.63 MeV. The new experimen- 
tal data have been used to determine the parameters of 
a locally optimized folding potential. In order to investi- 
gate the energy dependence of the potential parameters, 
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angular distributions of 89 Y(a,a) 89 Y elastic scattering at 
higher energies fl2l - fl6{ have also been analyzed. The vol- 
ume integrals of the local potential show a smooth vari- 
ation with energy over a broad energy range. However, 
the shape of the imaginary potential cannot be strictly 
fixed at low energies. 

In addition, the local potential is also used to study 
bound state parameters of the 93 Nb = 89 Y <S> a system. 
Excitation energies of excited states in 93 Nb and their 
decay properties can be described successfully with po- 
tential parameters very close to the scattering potential. 

The new experimental scattering data have also been 
used as a sensitive test for global parameterizations of 
the a-nucleus potential that have to be used in p pro- 
cess network calculations. Obviously, the best descrip- 
tion of the experimental data is obtained from the local 
potential where the parameters are fitted to reproduce 
the measured angular distributions. However, it is found 
that also the global parameterizations of [|| and pro- 
vide a good description for the measured angular distri- 
butions, and also reasonable agreement has been found 
for the angular distributions at slightly higher energies 
that are available from literature [12|, LL3| ■ For a deeper 
understanding of the differences between the locally opti- 
mized potential and the global a-nucleus potentials, the 
calculated elastic phase shifts Sl and scattered wave am- 



plitudes t]l are presented. 

Furthermore, the variation of elastic scattering angular 
distributions along the N = 50 isotonic chain has been 
analyzed. Here all global parameterizations failed to re- 
produce the amplitude and/or phase of the oscillations 
of the ratio of the Rutherford normalized cross sections. 
In order to advance our current understanding of the a- 
nucleus optical potential, further experimental scattering 
data with high precision are essential as well as improve- 
ments of the available global a-nucleus potentials. 
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